arXiv:1506.01554v2 [cond-mat.str-el] 27 Jun 2016 


Three magnons in an isotropic S = 1 ferromagnetic 
chain as an exactly solvable non-integrable system 


P.N. Bibikov 


Russian State Hydrometeorological University, Saint-Petersburg, Russia 

Abstract 

It is shown that a generalization of Bethe Ansatz based on an utilization of 
degenerative diserete-diffractive wave functions solves the three-magnon problem 


for the 5 = 1 isotropic ferromagnetic infinite chain. The four-magnon problem is 
briefly discussed. 


1 Introduction 

Quantum integrable models most usually are solved by various versions of Bethe Ansatz 


However solvability does not imply integrability |1]. In fact the later results in a 


non ergodic physical behavior, while the former gives the possibility to obtain the spectrum 
of physical states. An interplay between these conceptions may be well illustrated for 
models whose particles (elementary excitations) are flat waves excited from the ground 
state |0) by some creation operators Here the index j = parameterizes 

internal degrees of freedom such as polarization of triplons in spin ladders {d = 3) [5], or 
electrons in the t — J model (d = 2) [6]. Namely one- and two-particle states are 



( 1 ) 


n 



ni<n2 


1 



As vectors in C'^ 0 C'^ the two amplitudes and are related by the formula 

A 21 = S{ki,k2)Ai2^ (3) 

where the (P x cP matrix S{ki, k 2 ) is the two-particle scattering matrix. If it satishes the 
Yang-Baxter equation [7] 

Sl2{kh, kc)S23{ka, kc)Si2{ka, h) = S23{ka, h)Si2{ka, kc)S23{kb, /cj , (4) 

where (/ is the hnite-dimensional matrix unit) 

Si2ik,k) = Sik,k)®I, S23ik,k) = I^S{k,k), (5) 

then the wave function of a three-particle state 

\kuh,h)»»»= Y, ( 6 ) 

ni<n2<n3 

may be obtained in the Bethe form (we have omitted the polarization indices, and used 
the Levi-Civita tensor) 

3 

an„n2,n3{kl,k2,k3) = SabcAaUh, k2, (7) 

a,b,c=l 

which express both solvability and integrability. The former is evident, while the latter 
follows from the fact that according to (7) an initial set of incoming wave numbers does 
not change under the scattering. 

Equation (4) was hrst treated as an integrability condition for quantum gas with delta- 
function interaction [8] and then used as an integrability test for other models such as 
the above-mentioned spin ladder and t-J model mn- For a ID system equation (4) 
together with an absence of particle production are necessary conditions for factorization 
of the multi-particle scattering or equivalently for its reduction to a succession of space- 
time separated two-particle collisions [9]. For identical particles due to the energy and 
momentum conservation laws, such a collision reduces to an exchange of wave numbers 
between scattering particles and multiplication of amplitudes on an appropriate two- 
particle S'-matrix (see equation (3)). That is why an m-particle wave function is a sum 
of m\ exponential terms each of them corresponds to a permutation of wave numbers in 
a set {A:i,... ,km}- Since a distribution of wave numbers (momentums) does not alter 
under the collisions the system possesses a non ergodic behavior. 
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For a non-integrable quantum system the situation is drastically different. One- and 
two-particle states may be obtained as previously [HIE]. However a three-particle wave 
function should have a diffractive form 

„ 3 3 

am,n„n3 = / dhdk2dh6(^l[e^^‘ 

d 1=1 i=i 

3 

a,6,c=l 


corresponding to changes of an incoming triple of wave numbers and hence to an ergodic 
behavior. Unfortunately a substitution of (8) into (6) does not reduce the spectral problem 
to a rather simple form as it did in the case (7). 

In the present paper we suggest a discrete analog of the diffractive form (8), namely 


a 


ni,n2,n3 


M 3 


(9) 


Here M > 1 and 


1=1 


= const, 


3 

= const. 

1=1 


( 10 ) 


As will be shown here for a degenerative case M < oo the ansatz (9) results in essential 
simplifications. Namely for 

1< M <oo, (11) 

the system is still ergodic and non-integrable but the three-particle wave functions (9) 
may be derived exactly as in the usual version of Bethe Ansatz. 

Magnons (spin waves with AS = 1) in a 5 = 1 isotropic ferromagnetic chain (see the 
Hamiltonian (12)-(13) of the present paper) have no internal degrees of freedom. Hence 
the two-magnon S'-matrix is a scalar function and equation (4) is satisfied automatically. 
However the two-magnon scattering results in the creation of the quadruplon resonance 
[12] (the spin wave with AS* = 2). In the su(3)-symmetric point [13] the quadruplons 
are stable particles with the same energy as magnons. For a broken su(3)-symmetry the 
quadruplons become instable and decay into magnon pairs, at the same time creating as 
resonances in two-magnon collisions. Due to these processes a three-magnon scattering has 
the following channel. First of all two neighboring incoming magnons create a quadruplon 
resonance which then collide with the third incoming magnon. Under this collision the 
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resonance decays on two magnons. One of them goes to infinity while the other forms a 
new resonance with the third magnon. Finitely this new resonance again decays on two 
outgoing magnons. As it will be shown in Sect. 4 only this process can not be accounted 
by the wave functions of the form (7). However the M = 2 wave functions of the form 
(9) will be derived. 

2 Hamiltonian and one-magnon spectrum 

In the present paper we study the general model of isotropic S' = 1 ferromagnet with the 
Hamiltonian 

OO 

^ Hn,n+1, ( 12 ) 

n=—oo 

where [laiiij 

Hn,n+1 = i-{Sn- S.+l) + j(/ - (S„ • . (13) 

Here S„ is the standard triple of S' = 1 spin operators associated with n-th site. The 
constant terms proportional to the inhnite-dimensional matrix unit / are added only for 
the relation 

m = 0 , ( 14 ) 

where the state 

OO 

|0)= n + (15) 

n=—oo 

(the states \j)n with j = —1,0,1 form the standard S' = 1 triple associated with n-th 
site) will be treated as the pseudovacuum. It may be readily proved by an analysis of 
the spectrum of the 9x9 Hamiltonian density matrix related to Hn^n+i that for J < 1 
equation (15) gives the ground state of the model. For J > 1 it will be the ground state 
under the saturating magnetic held. 

Using the representation 

(s„. s „+0 = sf.s;„ +1 (s:s;+i + s;s:+,), (le) 

and relations 

S±|±1)„ = 0, S±|0)„ = x/2|±1)„, S±|Tl)n = \/2|0)„, Sl\j)n = 3\3)n. (17) 
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one readily gets [H] 


Hn,n+l \ i l)n ® | i l)n+l 0 

Hn,n+l\ ± l)n <8 |0)n+l = | ± 1)„ 0 |0)n+l “ |0)n <8 | ± l)n+l, 
i/n,n+l|0)n ® | ± l)n+l = |0)n ® | ± l)n+l “ | ± l)n ® |0)n+l, 

Hn,n+l\0)n ® |0)n+l = (1 — J) ^|0)n 0 |0)n+l ~ |l)n ® | — l)n+l ~ | ~ l)n ® |l)n+l^ , 

Hn,n+l\ ± l)n O I T l)n+l = (2 - J) | ± 1)„ 0 | ^ l)„+i + (J - 1)|0)„ 0 |0)n+i 

-J|Tl)n® |±l)n+l. (18) 


A one-magnon state is the flat wave 


whose energy 


l,A.) = ^e“”S.T|0>. 

n 


Emagnik^ 2(1 COs/c), 


may be readily obtained from (18). 


(19) 


( 20 ) 


3 The two-magnon scattering 


A two-magnon state should have the form na 


| 2 ) 


E 

ni<n2 


(^ni,n2^ni^n2 




( 21 ) 


According to (18) the corresponding Schrodinger equation on the combined wave function 
{ani,n 2 7 bn} has the form of the following system [H] 


^^ni,n2 —1,^2 ^ni+l,n2 ^ni,n2 —1 ^ni,n2+l -^^ni,n2? ^2 ^1 ^ I 5 

(3 l,n+l ^n,n-\-2 “1“ ('^ l)(^n “1“ ^n+l) EUn^n+li 

2(2 J^bn J(^bn—1 &n+l) “1“ l)(^n—l,n Ebn- 


( 22 ) 

(23) 


In the su(3) invariant point J = 1 the system (22)-(23) splits into two independent 
subsystems and a solution 

bn = 6 *"=", (24) 

with the same energy as ( 20 ) corresponds to the above-mentioned quadrupole wave [T^fT^ 
(quadruplon). Turning to the general case J 7 ^ 1 and suggesting the following ansatz 

an^n^iku k2) = bn{ki, k^) = (25) 
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we readily get from (22) the energy 


E{ki^ ^ 2 ) -^magn(^l) “1“ -^magn(^2)) (^6) 

and reduce (23) to 

1 + - (1 + J)e*^ 2 j 2I12 - - (1 + J)e^’^A A21 

+ (J- 1)(^1 + = 0, 

(J - 1) (^e-*^i + e*^^) ^12 - (J - 1) (e*^' + A 21 

+2^cos A:i + cos/c 2 — J(cos (/ci + k 2 ) + l)j-B = 0. (27) 

System (27) has the following solutions 

Ai 2 = A{ki,k2)i A 21 = A{k2iki), B = B{ki,k2), (28) 

where 

A{k,k) = e*^+e*(^+2fc)_^l^j)g2*fc_^^^3j^gi(fc+fc)^jj'3gifc^3gi(2fc+fc)_^_g2*(fc+fc)j^ ^29) 
and 

B{k, k) = {l- J) (^e*^ - e*^) (1 + e*(^+^)). (30) 

As we see the two-magnon problem may be solved in a non-diffractive way for all 
values of J. In other words it is unsensitive to non-integrability. Really equation (26) 
together with the condition ki + k 2 = k {k is the total momentum of the state) dehne the 
pair of wave numbers ki and /c 2 up to a permutation of them. 

4 The three-magnon problem 

A three-magnon state should have the form 
|3)=[ ^ (6W„(S-)^S- +(>g',„S-(S-)")]|0) (31) 

ni<n2<n3 ni<n2 

The corresponding Shrodinger equation on the combined wave function bnun2, ^^,112} 

splits on four groups of equations 

6flni ,112, ®ni-|-l,n2,n3 Cirii ,112+1,113 ®ni ,n2,n3-|-l ®ni —I,n2,n3 ,n2 —l,n3 (Ini ,112,113 — ! E0,^1,112,1131 

(32) 
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at n 2 — Til > 1, ns — n 2 > 1, 


(5 l,m,n+l 0,m—2,m,n 

~\~{J ~ 1 ) ^m-l,n ^rn,n ~ Eam-l,m,n^ 

(5-J) 

*^m,n,n+l *^m,n,n +2 ^m+l,n,n+l l,n+l l,n,n+l 


+(J -1) 


5 ( 2 ) ^ 5 ( 2 ) 


m,n ~ ^m,n+l 


= i?a 


'm,n,n+l 5 


ofs — 4- 

^ )^m,n ^m,n+l ^m,n—l ^ ^m—l,n ' ^m+ 1 , 


(1) ^ a(i) 


+ (J - 1) 


^m,m+l,n ~t“ Q-m—l,m,r 
r(2) ,(2) 


= 

^m,ni 


2(‘\ — — iP'^ — iP'^ — / P'^ + 

^V'-' ^ )^m,n ^m+l,n ^m—l,n ^m,n—l ‘ ^m,n+l 


+ {J - 1) 


= Eb^‘^'^ 

^m,n 5 


|^^m,n,nH-l H“ l,i 

at n — m > 1 , 

2(2 Q-n—l,n,n +2 Q-n— 2 ,n,n+l + ('^ 1) ^n—l,n+l “l“^n,n+l l,n“l“^n—l,n+l 

and 


(33) 


Edfi- 

(34) 


(4 ^i,i+2 J^nP,n+l ^i,n+l + ("^ l)fln-l,n,n+l — 

(4 ~ '^)^i-l,n ~ ^i-2,n ~ "^^n-l.n+l ~ ^n-l,n + ("^ ~ l)fln-l,n,n+l = (35) 

The following substitution 

3 

anun„n,{kl,k2,h) = eabcA{ka,k,)A{ka,K)A{h,k,)S^'^^^+^^^^^^^^^ 

a, 6 ,c=l 
1 ^ 

Pln,pip2p:,) = 2 ea6c5(A;a,fcfe)A(A;,,fce)dl(A:fe,A:e)e*'("“+"^)"^+"^"^], 

a, 6 ,c=l 

3 

CU(^1’ ^2, ^ 3 ) = 2 E ^-^cA{ka, h)A{ka, k,)B{h, (gg) 

a, 6 ,c=l 

where 2 l(A:, ^) and B{k, k) are given by equations (29) and (30), solves equations (32)-(34) 
giving the energy 

3 

E{ki, k2, ks) = ^ Emagn{kj). (37) 

i=i 

At the same time equation (35) turns into 


xV1(k^,hi, fcaje^^i+fe+ta)" = 0, j = 1,2, 


(38) 


l,n,nH-l 5 
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where 


1 

h, h) = - £abce^"^A{ka, k,) [(^E{ku k2, ks)-A + + j(l + e-'('^“+'='>))) 


a,b,c=l 


■B{ka, h)A{h, kc) + E^»A{ka, h)B{h, k,) + 2(1 - J)e-^^‘^A{ka, h)A{h, kc) 


x^^\k,, k2, ks) = ^ ^ Sabce-^’^'^Aika, k,) [(^E{k,, k^, k^) - 4 + + j(l + e‘('='>+'^=))) 

a,6,c=l 

■A{ka, h)B{h, k,) + e-^^'^Bika, h)A{h, fee) + 2(1 - jy^^A{ka, h)A{h, kc)] . (39) 


An evaluation of the sums in (39) with the use of the computer algebra system MAPLE 
gives 


xW(/ci, k2, /cs) = (1 - J^Mki, k2, ks) \(E{ki, k2, h) - S^e''^ - 1 + j(2 + 3e''= + 


J 2 ik 


k2, ks) = (1 - J^Mki, k2, ks) \e’^ - E{k^, k2, fca) + 5 - jfs + + e-'"^)! .(40) 


where k = ki + k 2 + k^ and 

3 

ip{k,, k2, ks) = (e'^i - e'^^) (e'^^ - e'^^) (e'^^ - ■ JJ(1 - e 


ikn 


(41) 


1=1 


As we see from (40) X^^^ki, k 2 , k^) = 0, {j = 1, 2) only in two integrable cases [131 [ 15 ] 

J = 1 and J = —1. However even in the general case the whole system (32)-(35) will be 
obviously satished for the following wave functions (j = 1,2) 

Q'ni,n2,n3 (/ci, /c2, ^3, ku k2, ks) = ip{ki , A)2, ^3)^ni,n2,?i-3 {ki, k2, ks) - V?(fci, k2, k3)an,,n2,n3{kl, k2, ks), 

^S,n2(^i> k2, ks, h, h, h) = h, h)bliln^{ki, k2, k^) - ip{ki, k2, k3)bl^l^^{ki, h, ^3),(42) 


where 


P[e* = n®""' E En^mCkl) = 


(43) 


1=1 


l=l 


l=l 


l=l 


5 Remark on the four-magnon problem 

An evaluation of four-magnon states for our model is a problem of the next level of 
complexity. In order to see this let us recall an evaluation of equation (42). First of 
all we take a triple {fci, ^2, ^3} and then construct the wave function (36) which satishes 
equations (32)-(34) but does not satishes equation (35). In order to solve the latter we 










add the term related to a new triple {/ci,/i: 2 ,/ca}. The resulting wave function has the 
form (42). 

For a four-magnon state 


1^) ^ ^ ^'i^l-,n2,ns^ri4^ni^n2^n^^n4 ^ ^ ^^ni,n2,n3 (^ni) ^n2^Ti3 

ni<n2<nz<n4 ni<n2<n^ 

+€’„,,»,s-(s-)"s- +(,g>„,,„,s-s-(s-)4 + c„„„,(s-)^(s;jd0)-(«) 


ni<n2 


the corresponding Schrodinger equation should again split on several systems of equations 
related to different processes in the four-magnon system. We shall study only one of them 
related to an extreme right magnon being separated from the others and hence do not 
interacting with them. 

As in the previous case we take at once ani,n 2 ,n 3 ,n 4 as a linear combination of 24 = 4! 
Bethe exponents related to different permutations of four wave numbers ki, /c 2 , ks and 
^ 4 . Let us hrst consider six of them proportional to From an account of the 

interaction between the triple of left magnons, it follows that these six terms should be 
added to another six ones proportional to the same exponent e*^"*”"* but with a new triple 
{ki, k 2 , ks} of the left magnons wave numbers. The same picture will be for all other three 
groups of exponents proportional to {j = 1, 2, 3). 

Now we may explain the cardinal difference between three- and four-magnon problems. 
In the former case for a given triple {ki, k 2 , k^} it is sufficient to add only a single triple 
{^ 1 , ^ 2 , ^ 3 }, however in the latter one for a given quartet {k^i\ k ^2 \ it is necessary 

to add at least four different new quartets 


fiLb) uy3)\ 


0> 


hii) _ jL{o) 




(45) 


related the same energy and total wave number. Each of the four induced quartets has 
the similar rights to the initial one and hence should be in the same correspondence with 
some other four quartets (one of them is the initial quartet). As a result the total set of 
quartets may be represented as a graph whose each vertex (related to its own quartet) is 
connected with four different other vertices. 

Let us study in detail a construction of the simplest example of such set of quartets. As 
it was explained above, hrst of all we take an initial quartet {k^i \ ■ ■ ■, kf^^} which induces 
the four new ones according to equation (45). Taking now the quartet {k^\ ..., k^^^} 
and applying the same argumentation as for {k ^\..., k^f^} we see that it also must be 
connected with four different quartets. One of them is already known: it is ..., kf^^}. 
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Hence, we should present the rest three quartets. The simplest way to do this is to use the 
quartets (45) with j = 2, 3,4 (otherwise we have to introduce new quartets {ki\ ..., k^^}, 
j = 5,6,7). Under this choice any pair {k[^\ ..., k^^^} and {ki\ ..., k^^^} {j = 2,3,4) 
should have a common wave number. According to equation (45) the wave numbers kj^^ 
(j = 1,..., 4) are already utilized. Hence without lost of generality we may put = k[^^ 
(j = 2, 3,4). Turning to the quartet {/c) ,... ,kl } we see that since it has been already 
connected with {kf’\ ..., kf^} and ..., k^^} we have to connect it only with two 

quartets. Again, the simplest way to do this is to use {ki \..., {j = 3,4). Since 
the wave numbers /c) with j = 1,2 are already utilized we may postulate (without any 
lost of generality) kf'^ = k^'^ {j = 3,4). Finitely we consider the quartet ..., kf ^} 

and connect it with {k ^\..., by the formula A:® = kl^\ As a result we have the 
system of five quartets (geometrically it may be represented as a graph with five vertices 
connected to each other) 

{k[^\...,k^^^}, i = 0,...,4, (46) 


and 10 relations; namely (45) and 

kl^^ = kf\ j,/ = l,...,4, (47) 

According to Eqs. (45) and (47) only 10 of the 20 wave numbers A:^^ (j = 0,.. .4, I = 
1,...4) for example A:^^ (0 < j < /, / = 1,...4) are independent. According to the 
energy and quasimomentum conservation laws they satisfy 10 independent equations 


+k^^^ +k^^ +k^^) 


= e 




= e 


= e 


(48) 


and 


4 4 

Si)\ _ 77 /;dl)N _ 77 


d2)^ 


^ Emagn^kj ) — E^agni^kj ) — -E'ma5n(^2 ) ^ -5'ynagn(^-y 

j=l j=l j=2 

3 4 

= ^ ^ Emagn{k^^) + Emagn{k^ ^) = ^ ^ Emagnik^i^) = E. (49) 

i=i i=i 

Correspondingly a general system of M quartets contains AM wave numbers. Since 
each of them is common to two different quartets, only a half of them (namely 2M) are 
independent. Equations 




( 50 ) 


1 = 1 


l=l 
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where k is the quasimomentum (total wave number) and E is the energy, give 2M con¬ 
ditions on these 2M wave numbers. Hence the existence of £nite-M four-magnon Bethe 
wave functions is in question even without an analysis of the pure four-magnon collisions. 


6 Summary and discussion 

In the present paper we have studied the three-magnon problem for a general isotropic S = 
1 ferromagnetic inhnite chain. Except the two integrable cases da [15] the corresponding 
wave functions can not be represented in the Bethe form (7) but only as a non-integrable 
modihcation (9), (42) (we call it the degenerative, discrete-diffractive form). Since the 
presented set of states is highly overloaded a complete description of the three-magnon 
scattering [16] may be obtained only after an extraction of a non-overloaded complete 
system of the three-magnon states. However, it is not clear how to represent such a system. 
In fact, for an integrable spin chain the three-magnon eigenstates may be parameterized 
by their energy E, quasimomentum k and the eigenvalue of an additional hrst integral. 
The latter belongs to an inhnite set of commuting hrst integrals which may be obtained 
by the standard procedure HI El El [IS]. The system studied in the present paper is 
however non-integrable. Probably there exist an operator Q which commutes both with 
the Hamiltonian and the shift operator. If its eigenvalue Q is independent from E and 
k, it may be used for a parametrization of the spectrum. Otherwise the parametrization 
procedure seems unclear and probably may be developed on the base of noncommutative 
geometry im. Nevertheless, it seems evident that a scattering of three incoming magnons 
with wave numbers {ki, k 2 , k^} should result in the creation of all possible outgoing three- 
magnon states with the same E and k. 

We also have shown that the corresponding four-magnon problem is much more diffi¬ 
cult. 

We suggest that the obtained result in its future development may be useful for a 
derivation of low-temperature expansions for thermodynamical quantities in the gapped 
regime [HI [H] [20] . 

Finitely we notice that although equation (4) on the S'-matrix has the same form as 
the equation on the so called i?-matrix (the Yang-Baxter equation in the braid group 
form [D El 0 [15]) the two subjects are not directly connected to each other. The S'- 
matrix characterizes a two-magnon scattering and its dimension (P x <P depends on the 
number of elementary excitations (that is d in our notations). From the other hand the 
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/^-matrix has the same dimension (P x (P as the Hamiltonian density matrix. Here d is 
the dimension of the Hilbert space associated with each site of the chain (for example 
d = A for spin ladders [5] and d = 3 for the t — J model 0). Usually d < d. Moreover 
as it was shown in the present paper the Yang-Baxter equation for the S'-matrix may be 
satisfied even in non-integrable cases when the i?-matrix formalism is irrelevant. 

The author is grateful to L. D. Faddeev for the helpful discussion and to H. Katsura 
for the useful comment. 
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